I> 



Revisiting the Phase Structure of the Polyakov-quark-meson Model in the presence of 

Vacuum Fermion Fluctuation. 

Uma Shankar Gupta* and Vivek Kumar Tiwari'l' 
Department of Physics, University of Allahabad, Allahabad 211002, India. 

(Dated: July 8, 2011) 

We have considered the contribution of fermionic vacuum loop in the effective potential of Polyakov 
loop extended Quark Meson Model (PQM) for the two quark flavour case and explored the phase 
structure and thermodynamics of the resulting PQMVT model (Polyakov Quark Meson Model with 
Vacuum Term) in detail at non zero as well as zero chemical potential. The temperature variations 
of order parameters and their derivatives have been calculated and the phase diagram together with 
the location of critical end point (CEP) has been obtained in /x, and T plane in both the models 
PQMVT and PQM. The PQMVT model analysis has been compared with the calculations in PQM 
model in order to bring out the effect of fermionic vacuum term on the physical observables. We 
notice that the critical end point (CEP) which is located near the temperature axis at (^ = 81.0, 
^— V \ T = 167 MeV) in the PQM model gets shifted close to the chemical potential axis at {^cep ~ 294.7, 

fvj ■ TcEP = 84.0 MeV) in the PQMVT model calculations of the phase diagram. Since it emerges from a 

background of second order transition in the chiral limit of massless quarks, the crossover occurring 
at /i = in PQMVT model for the realistic case of explicitly broken chiral symmetry, has been 
identified as quite a soft and smooth transition. We have presented and compared the results for 
temperature variations of thermodynamic observables at zero and different non-zero quark chemical 
potentials. It is noticed that the presence of fermionic vacuum term in the effective potential leads 
to a smoother and slower temperature variation of thermodynamic quantities. 

Q^' PACS numbers: 12.38.Aw, 11.30.Rd, 12.39.Fe, ll.lO.Wx 

I. INTRODUCTION 

J^ r Quantum Chromo-dynamics (QCD), the commonly accepted theory of strong interaction predicts that normal 

, ' hadronic matter undergoes a phase transition, where the individual hadrons dissolve into their constituents and 

fKf-N \ produce a collective form of matter knovifn as the Quark Gluon Plasma (QGP) under the extreme conditions of high 

temperature and/or density [1-4]. Relativistic heavy ion collision experiments at RHIC (BNL), LHC (CERN) and 

the future CBM experiments at the FAIR facility (GSI-Darmstadt) aim to create and study such a collective state 

i-^ I of matter. Study of the different aspects of this phase transition, is a tough and challenging task because Quantum 

y—^ ' Chromodynamics (QCD) becomes nonperturbative in the low energy limit. 

y—^ ,. It is well known that the basic QCD Lagrangian has the global SUL+R{Nf) x SUL~R{Nf) symmetry for Nf flavours 

L| ■ of massless quarks. The axial (A=L+R) part of this symmetry known as the chiral symmetry is spontaneously broken 

.^H , by the formation of a chiral condensate in the low energy hadronic vacuum of QCD and one gets (A^? — 1) massless 

^\. • Goldstone bosons according to the Goldstone's theorem. Since quarks are not massless in real life, chiral symmetry of 

5h I the QCD lagrangian gets explicitly broken and massless modes become pseudo-Goldstone bosons after acquiring mass. 

. . . 1 Nevertheless, the observed lightness of pions in nature suggests that we have an approximate chiral symmetry for 

QCD with two falvours of light u and d quarks. In the opposite limit of infinitely heavy quarks, QCD becomes a pure 

SU{Nc) gauge theory which remains invariant under the global Z{Nc) center symmetry of the color gauge group. The 

Center symmetry which is a symmetry of hadronic vacuum, gets spontaneously broken in the high temperature/density 

regime of QGP. The expectation value of the Wilson line (Polyakov loop) is related to the free energy of a static color 

charge. It vanishes in the confining phase as the quark has infinite free energy and becomes finite in the deconfined 

phase. Hence the Polyakov loop serves as the order parameter of the confinement-deconfinement phase transition [5] . 

Even though the center symmetry is always broken with the inclusion of dynamical quarks in the system, one can 

regard the Polyakov loop as an approximate order parameter because it is a good indicator of a rapid crossover in the 

confinement-deconfinement transition [6, 7]. 

The first principle lattice QCD Monte Carlo simulations (see e.g. [8-17]) give us important information and insights 
regarding various aspects of the QGP transition, like the restoration of chiral symmetry in QCD, order of the 
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confincment-dcconfincment phase transition, richness of the QCD phase structure and mapping of the phase diagram. 
Unfortunately progress in lattice QCD calculations has got severely hampered due to the QCD action becoming 
complex on account of the fermion sign problem [8] when baryon density/chemical potential is non zero. Though 
several methods have been developed to circumvent the sign problem at small baryon chemical potentials, a general 
solution to the sign problem for all chemical potentials is yet to be devised. Further since lattice calculations are 
technically involved and various issues are not conclusively settled within the lattice community, one resorts to the 
calculations within the ambit of phenomenological models developed in terms of effective degrees of freedom. These 
models serve to complement the lattice simulations and give much needed insight about the regions of phase diagram 
inaccessible to lattice simulations. 

In recent years, effective chiral models, having the pattern of chiral symmetry breaking as that of QCD like 
the linear sigma models(LSM) [18-24], the quark-meson (QM) models(see e.g. [25-35]), Nambu-Jona-Lasinio (NJL) 
type models [25, 36-39], have led to the investigation of the properties and structure of chiral symmetry restoring 
phase transition at sufficiently high temperature and density. Further these models were extended to incorporate 
the features of confinement-deconfinement transition where chiral condensate and Polyakov loop got simultaneously 
coupled to the quark degrees of freedom. Thus Polyakov loop augmented PNJL models [40-56] ,PLSM models and 
PQM models [57-65] have facilitated the investigation of the full QCD thermodynamics and phase structure at zero 
and finite quark chemical potential and it has been shown that bulk thermodynamics of the effective models agrees 
well with the lattice QCD data. 

In most of the QM/PQM model calculations, the fermion vacuum contributions to the free energy is frequently 
neglected[25, 26, 31, 33, 38] because here, the spontaneous breaking of chiral symmetry is generated by the mesonic 
potential itself. While in the NJL/PNJL model investigations, fermion vacuum term leads to the dynamical breaking 
of the chiral symmetry, hence it gets explicitly included up to a momentum cutoff A. Very recently, it has been 
shown by Skokov ct. al. in Ref. [64] that in a mean field approximation, where the fermion vacuum contribution 
to the free energy is neglected, the order of the phase transition for two flavour QM model in the massless chiral 
limit becomes first order at zero baryon chemical potential. They have further shown that the quark-meson model, 
with appropriately renormalized fermionic vacuum fluctuations in the thermodynamic potential, becomes an effective 
QCD-like model because now it can reproduce the second order chiral phase transition at /i = as expected from 
the universality arguments [66] for the two massless flavours of QCD. It has also been shown that in the presence 
of an external magnetic fleld, the structure of the phase diagram in the PQM model is considerably affected by the 
fermionic vacuum contribution [67]. In this paper, we will investigate the effect of fermionic vacuum fluctuations on 
the phase structure and thermodynamics of PQM/QM models in detail at non zero as well as zero chemical potential. 
In order to bring out the effect of fermionic vacuum term on the physical observables, we will compare the results of 
our calculation with the corresponding PQM model calculations without vacuum term. 

The arrangement of this paper is as follows. In Sec. II, we have given the formulation of PQM model for the two 
quark flavour. The Polyakov loop potential and the thermodynamic grand potential has been given in subsection 

II A. After giving a brief description of the appropriate renormalization of fermionic vacuum loop contribution, the 
subsection II B describes how the new model parameters are obtained in vacuum when renormalized vacuum term 
is added to the effective potential. The section III investigates the effect of fermionic vacuum term on the phase 
structure and thermodynamics. The subsection III A explores how, the temperature variation of order parameters 
and their derivatives at different chemical potentials, structure of the phase diagram in the /i and T plane and the 
location of critical end point, gets affected in the presence of vacuum term. The effect on the temperature variation 
of thermodynamic observables namely pressure, entropy, energy density and interaction measure has been discussed 
in the subsection III B while the discussion of specific heat, speed of sound and 2^4 has been presented in subsection 

III C and finally the subsection III D describes the results for quark number density and quark number susceptibility. 
Summary together with the conclusion has been presented in Sec. IV. The first and second partial derivatives of Uiog 
and rjqq with respect to temperature and chemical potential has been evaluated in appendix A. 

II. MODEL FORMULATION 

We will be working in the two flavor quark meson linear sigma model which has been combined with the Polyakov 
loop potential [57]In this model, quarks coming in two flavor are coupled to the SUl{'2) x SUii{2) symmetric four 
mesonic fields a and tt together with spatially constant temporal gauge field represented by Polyakov loop potential. 
Polyakov loop field $(af) is defined as the thermal expectation value of color trace of Wilson loop in temporal direction 

$ = — Tr^L, $* = — Tr^Lt (1) 



where L(x) is a matrix in the fundamental representation of the S'f/c(3) color gauge group. 
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Here V is path ordering, ^o is the temporal component of Euclidean vector field and (3 ^ T ^ [5] . 

The model Lagrangian is written in terms of quarks, mesons, couplings and Polyakov loop potential U ($, $*, T). 

Cpqm^Cqm-U{^,^\T) (3) 

where the Lagrangian in quark meson linear sigma model 

^QM = Qf [il^D^, - g{a + i'^^f ■ tt)] qj + Cm (4) 

The coupling of quarks with the uniform temporal background gauge field is effected by the following replacement 
D^i = dfj, — iAfj, and A^ = <5^o^o (Polyakov gauge), where A^ = gs^JiA'^/2. gg is the SUc{3) gauge coupling. Xa are 
Gell-Mann matrices in the color space, a runs from 1 • • • 8. qj — (u, d) denotes the quarks coming in two flavors and 
three colors, g is the flavor blind Yukawa coupling that couples the two flavor of quarks with four mesons; one scalar 
(ct, J^ = 0+) and three pseudoscalars (tt, J^ =0^). 

The quarks have no intrinsic mass but become massive after spontaneous chiral symmetry breaking because of 
nonvanishing vacuum expectation value of the chiral condensate. The mesonic part of the Lagrangian has the following 
form 

Cm - \{d,af + ]^{d,iTf-U{a,n) (5) 

The pure mesonic potential is given by the expression 

t/(a,^) = ^((j2 + 7r2_^,2^2_^^^ ^g^ 

Here A is quartic coupling of the mesonic flelds, v is the vacuum expectation value of scalar field when chiral symmetry 
is explicitly broken and h =/^m^ . 

A. Polyakov loop potential and thermodynamic grand potential 

The effective potential U ($, $*,r) is constructed such that it reproduces thermodynamics of pure glue theory on 
the lattice for temperatures upto about twice the deconfincment phase transition temperature. In this work, we are 
using the logarithmic form of Polyakov loop effective potential Ref. [41]. The results produced by this potential are 
known to be fitted well to lattice results. 

+4($*3 + $3)-3($*$)2] (7) 

where the temperature dependent coefficients are as follow 

a(r) = ao + «,(!) +«,(!)' 6(T) = 63(f 

The critical temperature for deconfinement phase transition Tg = 270 MeV is fixed for pure gauge Yang Mills 
theory. In the presence of dynamical quarks Tq is directly linked to the mass-scale Aqcd , the parameter which has 
a flavor and chemical potential dependence in full dynamical QCD and Tq — > To(7V/,/Lt) [57, 65]. For our numerical 
calculations in this paper, we have taken a flxed Tq ~ 208 for two flavours of quarks. 

The parameters of Eq.(7) are 

ao = 3.51, ai = -2.47 , 

02 = 15.2 , &3 = -1.75 



In the mean-field approximation, the thermodynamic grand potential for the PQM model is given as [57] 

C!mf(T, ^; a, $,$*)== U{T- $,$*) + U{a) + r!,,-(T, /^; a, $, $*). 

Here, we have written the vacuum expectation values (cr) = a and (tt) — 
The quark/ ant iquark contribution in the presence of Polyakov loop reads 



(8) 



C!,,(r, /i; a, $, $*) = n^f + OT = -27V/ j -^ {N,E,e{K^ - f) + r[ln.g+ + In.g-] } 



(9) 



The first term of the Eq. (9) denotes the fermion vacuum contribution, regularized by the ultraviolet cutoff A. In 
the second term gj" and g~ have been defined after taking trace over color space. 
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1 + 3$e" 



-K/T + 3$*e-2^'^/'^ + e-^^i'^ 



(10) 



1 + 3$*e-^."/^ + 3$e-2S,-/T + g-3£-/T 



(11) 



Here we use the notation E^ ~ Eq^ ^ and Eq is the single particle energy of quark/antiquark. 



Er, 



p^ + rUq^ 



(12) 



where the constituent quark mass rUq = ga is a function of chiral condensate. In vacuum cr(0, 0) = ctq = A = 93.0MeV 

B. The renormalized vacuum term and model parameters 

The fermion vacuum loop contribution can be obtained by appropriately renormalizing the first term of Eq. (9) 
using the dimensional regularization scheme, as done in Ref.[64]. A brief description of essential steps is given below. 
Fermion vacuum term is just the one-loop zero temperature effective potential at lowest order [68] 






-2NfN, I S^E 



(27r) 



3^q 



= -2NfN, 



/(^l^^(^" + ^^')+^' 



(13) 



the infinite constant K is independent of the fermion mass, hence it is dropped. 

The dimensional regularization of Eq. (13) near three dimensions, d = 3 — 2e leads to the potential up to zeroth 
order in e as given by 



«■? 167r2 Me 2 



-^ + ^-- + "°(27S7) 



(14) 



here M denotes the arbitrary renormalization scale. 

The addition of a counter term SC in the Lagrangian of the QM or PQM model 
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gives the renormalized fermion vacuum loop contribution as 
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(15) 



(16) 



Now the first term of Eq. (9) which is vacuum contribution will be replaced by the appropriately renormalized 
fermion vacuum loop contribution as given in Eq. (16). 



The relevant part of the effective potential in Eq. (8) which will fix the value of the parameters A and v in the 
vacuum at T = and /i = is the purely a dependent mesonic potential U{a) plus the renormalized vacuum term 
given by Eq. (16). 

^cNf 4 4 (g(T\ Aw2 A 

The first derivative of ^{a) with respect to a at a ~ f^^ in the vacuum is put to zero 

af7MF(0,0;a,$,$*) dn{<j) 



nia) = 17- + Uia) = -^/.^ In (!^) - ^.^ + >^ - ha, (17) 



da da 

The second derivative of n{a) with respect to a at a = fj^ in the vacuum gives the mass of a 

2 _ d^nMF{0,0;U,^,<^*) _ dMa) 
"'" da^ da^ 

Solving the equations (18) and (19), we obtain 



(18) 



N,Nf 4 



3 + 41n,^^ 



(19) 



(20) 



and 



.,2_r^.,2, , NrNf , . 



Xv^ = {\v%+--^g^fi (21) 



where A^ and {Xv )a are the values of the parameters in the pure sigma model 



A. - 2/^ (22) 

(Az;^). = !!i_^ (23) 

It is evident from the equations (20) and (21) that the value of the parameters A and v^ have a logarithmic 
dependence on the arbitrary rcnormalization scale M. However, when we put the value of A and Au^ in Eq.(17), the 
M dependence cancels out neatly after the rearrangement of terms. Finally we obtain 

8^2 y \fj 2 4 

Here, wc define A^ and A^v^ as the values of the parameters after proper accounting of the renormalized fermion 
vacuum contribution. 

A. = A. + ^/ (25) 



n{a) = -1-^gV^ In ^ - ^(r' + -^a^ - ha, (24) 



and 



,.2 .^..2^ , NcNf , ^2 



Kvt = i\v% + ^g^ /^ (26) 

Now the thermodynamic grand potential for the PQM model in the presence of appropriately renormalized fermionic 
vacuum contribution (PQMVT model) will be written as 

QuFiT, fi; a, $, $*) - W(r; $, $*) + n{a) + n^^{T, fi; a, $, $*). (27) 

Thus in the PQMVT model, One can get the chiral condensate a, and the Polyakov loop expectation values $, $* 
by searching the global minima of the grand potential in Eq.(27) for a given value of temperature T and chemical 
potential 

We will take the values ttItt = 138 MeV, rricr = 500 MeV, and /^ = 93 MeV in our numerical computation. The 
constituent quark mass in vacuum ?7i° = 335 MeV fixes the value of Yukawa coupling g — 3.3. 



III. EFFECT OF THE VACUUM TERM ON THE PHASE STRUCTURE AND THERMODYNAMICS 

We are presenting the results of our calculation for studying the temperature variation of the order parameters a^ 
$, $*, their temperature derivatives and various thermodynamic observables at zero and non zero quark chemical 
potentials in the presence of the renormalized fermionic vacuum term in the effective potential of the PQM model. 
These results have been termed as PQMVT model calculations and we have investigated the interplay of chiral 
symmetry restoration and confinement-deconfinement transition in the influence of fermionic vacuum term. The 
phase diagram together with the location of critical end point (CEP) has been obtained in /^, and T plane for both 
the cases with and without fermionic vacuum contribution in the effective potential. In order to have a comparison, we 
have also shown the temperature variations of order parameters and their derivatives in the PQM model calculation 
with the same parameter set. The temperature variations of thermodynamic observables namely pressure, energy 
density and entropy density at three different chemical potentials (zero, [icEP and \x > jjlcep) have been shown in 
PQMVT model calculations. In order to study the effect of fermionic vacuum term at zero chemical potential, the 
temperature variation of the interaction measure, speed of sound, p/e ratio and specific heat, has been calculated in 
PQMVT model and QMVT model (Quark Meson model with vacuum term) and these results have been compared 
with the corresponding results in the PQM and QM model calculations. Finally we will be presenting the results of 
the temperature variation of baryon number density and quark number susceptibility at different chemical potentials 
in PQMVT model calculation. 

A. Phase structure 

The solutions of the coupled gap equations, Eq.(28) determine the nature of chiral and deconfinement phase 
transition through the temperature and chemical potential dependence of chiral condensate (T,the expectation value 
of the Polyakov loop $ and $*. Fig. 1(a) shows the temperature variation of the chiral condensate a normalized with 
the vacuum value on the left while the right end of the plot shows the Polyakov loop $ and $* temperature variation 
for the PQMVT model calculations, the corresponding temperature variation of the chiral and Polyakov loop order 
parameters in PQM model calculations has been shown in Fig. 2(a). In Fig. 1(a), the continuous dots, thin dash and 
thin solid lines represent the variation of — on the left and $ on the right at /i = 0, 294.7 (CEP) and 300 MeV 
respectively. Thick dash and thick solid lines in the right end of the plot represent the $* variations at ^ = 294.7 
and 300 MeV respectively. Fig. 1(b), 1(c) and 1(d), show the temperature derivatives of cr, $ and $* fields as a 
function of temperature respectively at three different chemical potentials ^ = ,100 and 294.73 MeV in PQMVT 
model calculations while the temperature variations of the same derivatives in the PQM model at /i = McV has 
been shown in Fig. 2(b). The characteristic temperatures (pseudocritical temperature) for the chiral transition T^ and 
the confinement-deconfinement transition T*, are defined by the peak positions (infiection point) in the temperature 
derivatives of a and $, $* fields. 

The chiral crossover transition for the realistic case of explicitly broken chiral symmetry, becomes quite soft and 
smooth at /i = because the corresponding chiral phase transition for massless quarks turns second order in the chiral 
limit after having a proper accounting of the fermionic vacuum contribution in the PQMVT model. The smoothness 
of crossover at ^ = is evident from the temperature variation of the chiral order parameter in Fig. 1(a) while the 
Polyakov loop order parameter variation at the same chemical potential, is sharp in comparison. The chiral crossover 
at /i = becomes less smooth as we increase the chemical potential. We find a large range (/i = at T^ = 186.5 MeV 
to ^cEP = 294.7 McV at T^ = 84.0 McV) in the values of chemical potential that makes the temperature variation 
of chiral order parameter, sharp and sharper such that eventually the crossover turns into a second order transition at 
CEP. The narrow width of the coincident variation of $ and $* temperature derivative at zero chemical potential in 
Fig. 1(b), signifies a sharp crossover for the confinement-deconfinement transition at T* = 169.0 MeV The a derivative 
shows a broad double peak structure at ^ = similar to the findings of N JL model calculation in Ref [48] , we have 
identified the chiral crossover temperature T^ = 186.5 MeV as the second peak position at higher pseudocritical 
temperature in Fig. 1(b). The first peak in the a derivative is driven by the sharp peak of the Polyakov loop variation. 
As the chemical potential is increased, the variation of Polyakov loop $ derivative becomes smoother and broader 
with increasing width, while the a derivative variation shows a decreasing width and double peak structure starts 
getting smeared after ^ — 100 MeV as shown in Fig. 1(c). For the chiral crossover transition in the chemical potential 
range /^ — 100 to 160 MeV, the identification of pseudocritical temperature T^ becomes ambiguous with an ambiguity 
of about 5 MeV. For /i > 160 MeV in the PQMVT model, double peak structure disappears from the temperature 
variation of the chiral order parameter temperature derivative as shown in Fig. 1(d) and its width decreases becoming 
narrow, narrower and narrowmost till the CEP at /i = 294.7 MeV and T = 84.0 MeV is reached where the chiral 
transition turns second order. 

For the realistic case of explicitly broken symmetry, the temperature variation of chiral order parameter at /Lt = 0, 
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FIG. 1: Temperature variations in the PMQVT model, (a) The continuous dots , dash and sohd lines represent the variation 
of ^ on the left end and $ on the right end of the plot at /x = 0, 294.73 (CEP) and 300 MeV respectively. Thick dash and 
thick solid lines in the right end of the plot represent the $* variations at /i = 294.7 and 300 MeV respectively, (b), (c) and 
(d), show the temperature derivatives of a, $ and $* fields as a function of temperature respectively at three different chemical 
potentials fi= 0, 100 and 200 MeV. 



turns out to be quite sharp and rapid in Fig. 2(a) in comparison to the corresponding PQMVT model variation in 
Fig. 1(a) because the chiral transition in the massless quark limit, is first order in the PQM model. Further the chiral 
transition remains a crossover in quite a small range from // = at T^ = 171.5 MeV to /i = 81 MeV at T^ = 167 
MeV in the PQM model results of Fig. 2(a). Since the chiral crossover is sharper than the confinement-deconfinement 
crossover in the PQM model calculations, the single peak of the a field temperature derivative in Fig. 2(b) at /i = 
is narrower and a lot higher than the peak in the variation of temperature derivatives of $ and $*. We have scaled 
the variation of $ and $* temperature derivatives in Fig. 2(b) by a multiple of 5 which shows a very small double 
peak kind of structure. We consider the chiral and confinement-deconfinement crossovers nearly coincident at /i = 0, 
T^ = 171.5 and we get exact coincidence as we move towards the CEP (T = 167.0 MeV and fj, = 81.0 MeV) of the 
model where on account of the transition turning second order, we get highest and narrowmost peak. 

In order to probe the issue of double peak structures emerging in Fig. 1(b), 1(c), the temperature derivatives of 
a, $ and $* fields have been evaluated as function of temperature by taking the Polynomial form [40] for Polyakov 
loop potential instead of the logarithmic form in PQMVT model. AT fi ^ 0, none of the field derivatives have 
double peak structure. It starts appearing at /i ~ 200 MeV separately in $ and $* derivatives as shown in Fig. 3(a) 
and we find robust noncoincident second peaks respectively at T* = 176 MeV, T* =156 MeV for fi = 280 MeV 
as shown in Fig. 3(b). The first and highest peak in a field temperature derivative gives the location of CEP at 
TcEP = 78.2 MeV and hcep = 293 MeV. Though we have not evaluated the whole phase diagram for this case, 
we find that the confinement-deconfinement transition line (obtained from nearly coincident peaks in $ and $*) 
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FIG. 2: (a) The continuous dots, dash and sohd hues in the left half of the figure represent the variation of -2- in the PQM 
model at /i = 0, ^ = 81 and /i = 130 MeV respectively. In the right end of the plot, continuous dots represent coincident 
variation of $ and <!>* at /i = while thick and thin dash lines represent the $* and $ variations at /i = 81 MeV respectively, 
(b) shows the temperature derivatives of a, $ and $* fields as a function of temperature at /i = in the PQM model. 



lies below the chiral crossover transition line (T* < T^) in the chemical potential range /x = to /i = 200 MeV 
and confinement-deconfinement crossover transitions for $* and $ fields for /x > 200 MeV separate in constituting 
different lines which get located above the chiral crossover phase boundary from jjl = 200 to /x = hcep = 293 MeV 
. These findings being similar to the results of Ref.[38, 65] support the quarkyonic phase[70] like scenario, having a 
region of confinement with restored chiral symmetry. It has been argued in Ref.[38] that transition temperature in 
the chiral sector decreases as the chemical potential is increased while the remnant of the deconfinement transition 
remain uneffected by the value of chemical potential. We notice that this explanation does not work in a calucation 
with logarithmic ansatz for Polyakov loop potential. Thus the meaning of the peaks in $, $* and a field temperature 
derivatives and their relation to the nature of crossover transition in confinement-deconfinement and chiral sector 
is debatable. Similar to the findings of NJL model calculation in Rcf[48], our calculations in PQMVT model 
with phenomenologically improved logarithmic Polyakov loop potential (which describes the gluon dynamics more 
appropriately with a better correlation to the effect of dynamical quarks) , show double peak in a field temperature 
derivative while the Polyakov field $ derivative shows a single peak in a chemical potential range /x = to /x = 200 
MeV as discussed in the second paragraph. Here we find support for the standard scenario [57, 65] where chiral 
symmetry restoration occurs at a higher temperature than the deconfinement transition. Thus the quarkyonic phase 
scenario is ruled out in our PQMVT model calculation with constant Tg when we have taken a logarithmic form for 
the Polyakov loop potential while the calculation with polynomial Polyakov loop potential supports its occurrence in 
certain range of fi and T. 

In Fig. 4, we have obtained the phase diagram in our calculation with logarithmic Polyakov loop potential and 
located the critical end point (CEP) in the PQMVT as well as PQM model calculations for rricr = 500 MeV. The 
structure of the phase diagram is very sensitive to the chosen value of the sigma meson mass. For the value ma- = 600 
MeV, the transition becomes a crossover in the entire /z and T plane for the PQMVT model calculation. We have 
shown the chiral crossover transition by a dash line which starts from T^ = 186.5 MeV at y^t = axis and ends at 
CEP; TcEP = 84 MeV and hcep = 294.7 MeV in PQMVT model. Due to the smearing of double peak structure in 
the temperature derivative of chiral order parameter in the range fi ~ 100 to 160 MeV, the chiral crossover transition 
temperature T^ is identified with an ambiguity of about 5 MeV. We get a unique T/ for /x > 160 MeV in the phase 
diagram because of a single peak structure which gets narrow and narrower for higher chemical potentials till we reach 
the CEP. The dash dotted line which starts at T* = 169 MeV and ends at CEP of the PQMVT model, signifies the 
confinement- deconfinement crossover transition. The chiral and confinement-deconfinement crossover transition lines 
merge at fi — 250 and T^ = T* = 132. The thin solid line for /i > i^cep represents the first order phase transition 
corresponding to the jump in all the order parameters at the same critical temperature. The chiral crossover transition 
line lies above the crossover line for the confinement-deconfinement transition. Thus our results of the PQMVT model 
calculation are in tune with the standard scenario [65] where chiral symmetry restoration occurs at a higher critical 
temperature T^ = 186.5 MeV than the confinement- deconfinement transition temperature Tf = 169 MeV at /.i = 
axis. Further the crossover transition temperature at /z = compare well with the lattice [12, 15] results and QCD 
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FIG. 3: Temperature variation of order parameter derivatives with polynomial Polyakov loop potential in PQMVT model 
calculation. 
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based computations [65, 69] in two flavour model. The chiral and confinement-dcconfincmcnt crossover transition 
lines are coincident (as shown by the continuous dots) and start from T^ = Tf = 171.5 MeV at // = MeV to end at 
the CEP of the PQM model. The first order transition for /i > hcep in the PQM model, has been shown by the thick 
solid line. The CEP of the PQM model gets located near the temperature axis at ficEP = 81 MeV and Tqep — 167 
MeV because the chiral crossover at ^ = 0, having the background of a first order phase transition in the chiral limit, 
is rapid and sharp and soon it gets converted to a first order phase transition as we increase the chemical potential. 
While the critical end point (CEP) gets shifted close to the chemical potential axis in PQMVT model because the 
chiral crossover transition at ^ = is quite soft and smooth as it emerges from a phase transition which turns second 
order in the chiral limit due to the effect of renormalized fermionic vacuum contribution in the effective potential and 
further it remains a crossover for large values of the chemical potential. 

B. Thermodynamic ObservablesrPressure, Entropy and Energy Density 

The negative of grand potential gives the thermodynamic pressure 



PiT./J-) = -riMF(T,/z) 



(29) 
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Thermodynamic pressure divided by the QCD Stcfan-Boltzmann (SB) hmit has been shown for three chemical 
potentials ^ = 0,294.7 (CEP) and 300 MeV in Fig.5 for PQMVT model. It has been normalized to vanish at 
T = n = 0. We have shown the the pressure calculated in PQM model also for comparison at ^ = 0. For Nf massless 
quarks and N^ — 1 massless gluons in the deconfined phase, the QCD pressure in the SB limit is given by 



PSB 

jn4 



{N^-l)-+N,Nf 



180 "^6 



6 \T/ 127r2 \tJ 



(30) 



The fermionic vacuum contribution makes the pressure variation in PQMVT model smooth at fi — and this 
curve (thin solid line) lies slightly below the curve (line with continuous dots) obtained in PQM model. The pressure 
variations at ^cep = 294.7 and n = 300 MeV of PQMVT model are represented by the thick solid and dash line 
respectively.The pressure increases near the chiral transition due to the melting of the constituent quark masses and 
saturates at about eighty percent of the SB limit. 

The entropy density is defined as negative of the temperature derivative of the grand potential. 



dn 



MF 



dT 



(31) 



The implicit variation of cr, $ and $* fields with respect to temperature has been accounted for, in the temperature 
derivative of ^{(J)^ U\og and fi^g as evaluated in the appendix A. The temperature variation of entropy density 
normalized by its QCD SB limit has been shown in Fig. 6. It is continuous for crossover transition and attains about 
40-45 percent of its SB value at psedocritical transition temperature. Again due to the fermionic vacuum fluctuations, 
the entropy density variation (thin solid line) at /x = turns out to be a smoother function of temperature in 
PQMVT model when it is compared with corresponding curve (line with continuous dots) of PQM model calculation. 
At /icEP — 294.7 MeV, the entropy density curve (thick solid line) shows a steep rise at Tqep = 84. MeV in PQMVT 
model, then it takes a bend to reach its saturation. The PQM model entropy density curve (dah dotted line) at 
/i = 294.7 MeV shows a large jump because chiral transition is strong first order at this chemical potential. The first 
order chiral transition of PQMVT model at // = 300.0 MeV, generates another jump in the entropy density curve 
(line with dash), though this jump is smaller than the first order jump seen in PQM model entropy curve at a lower 
chemical potential /x = 294.7 MeV. 

The energy density in the presence of chemical potential is given as 



e = — p + Ts + ^n 



(32) 



where n is the number density.The temperature variation of energy density normalized by its QCD SB limit value has 
been shown in Fig.7 for ^i = 0, 280 , 294.7 (CEP) and 300 MeV in PQMVT model. The energy density variation (thin 
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FIG. 6: Entropy variation with respect to temperature. 
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FIG. 7: Variation of energy density with respect to temperature. 



solid line) at fi — similar to entropy density variation, is smoother in comparison to the corresponding variation in 
PQM model calculation (line with continuous dots), this again is due to the influence of fermionic vacuum fluctuations. 
Similar to the entropy density variation at ^cep = 294.7 MeV, the energy density (thick solid line) also shows a very 
steep and large rise at Tcep = 84.0 MeV, then it curves to attain the saturation. At /i = 300.0 McV, we get a large 
jump in the energy density curve (dash-dotted line) which of course is a signature of the first order chiral transition. 
Since quark degrees of freedom get liberated and become light, the energy density registers a rapid increase near the 
crossover/phase transition point and reaches almost to the value of SB limit. 

The trace anomaly of energy momentum tensor is also known as interaction measure. The temperature variation 
of the interaction measure A = {E - 3p)/T^ has been shown in Fig.8 at /^ = MeV in PQMVT, QMVT and PQM, 
QM model calculations. The QM model variation of the interaction measure (line with continuous dots) shows a 
sharp and narrow peak near the pseduocritical transition temperature which becomes very broad and smooth in the 
corresponding variation ( thick dash line) of QMVT model calculation due to the effect of inclusion of fermion vacuum 
term contribution in the effective potential of QM model. The peak of interaction measure temperature variation 
(solid line) in PQMVT model shifts to a slightly higher temperature value in comparison to the corresponding peak 
in the variation (line with thin dash) of PQM model calculation. 
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C. Specific heat Cv and Speed of sound Cs 
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FIG. 9: Specific heat variation with respect to temperature. 



The expression of specific lieat at constant volume is given by 






= -T 



d^n 



MF 



QT2 



(33) 



The second partial temperature derivatives of a, $ and $* fields contribute in the double derivatives of fl{a), 
hl\og and riJq with respect to temperature as given in the appendix A. Fig. 9 shows the temperature variation of the 
specific heat Cy normahzed by T^ in QM, QMVT and in PQM, PQMVT model calculations at /i = 0. The specific 
heat variation while growing with the temperature, peaks at the crossover transition temperature and then saturates 
at the corresponding SB limit at the higher temperature. The QM model specific heat variation shows a large and 
sharp peak which becomes quite smooth and broad in the corresponding variation of QMVT model calculation due 
to the presence of fermionic vacuum term, further the peak position gets shifted to a higher transition temperature. 
The qualitative difference of structures in the curves of QM and QMVT model gets reduced due to the influence 
of Polyakov loop potential and we notice that the PQM model specific heat variation has quite a high and sharp 
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peak which becomes small and a little less sharp in the PQMVT model variation and the peaks occur at the same 
transition temperature. Further, we remark that the peak positions of the temperature variation of order parameter 
derivatives in Fig. 1(b) and Fig. 2(b) give different transition temperatures for chiral crossover in PQM and PQMVT 
model calculations while for confinement-deconfinement crossover, the transition temperature is almost same in both 
the models. 

The speed of sound is an important quantity for hydrodynamical investigations of relativistic heavy-ion collisions. 
It is given by 



C 



dp 


dp 




/de 


de 


5" 9T 


vl 


dT 



Cv ' 



(34) 



The equation of state parameter p(T)/e{T) also represents the information contained in trace anomaly. The velocity 
of sound Cg and the equation of state parameter p{T)/e(T) ratio has been shown as a function of temperature in 
QM, QMVT and PQM, PQMVT model calculations at ^ = in Fig. 10. Thick lines denote the result for the sound 
velocity C^ and thin lines show the variation of p{T) / e(T) ratio. The presence of fermion vacuum term in QMVT 
model leads to a very smooth temperature variation for C^ (line with thick long dash) and p{T)/e{T) ratio (line with 
thin long dash). The C^ temperature variation (line with thick continuous dots) in the QM model calculation, shows 
a very sharp drop followed by a rapid rise while the EOS parameter p{T) / e(T) ratio (line with thin continuous dots) 
shows a cusp at crossover transition temperature. The temperature variation of C^ (thick solid line) and p{T) / e{T) 
ratio (thin solid line) in the PQMVT model turns out to be smoother than the corresponding variation of C^ (line 
with thick, short and dark dash) and p{T)/e{T) ratio (thin dash line) in the PQM model calculation. At higher 
temperatures C^ and p{T)/e{T) ratio approach the ideal gas value 1/3 in all the cases of model calculation. In PQM 
and PQMVT models, the value of C^ almost matches with the p{T)/e{T) ratio for lower temperatures and C^ value 
is always larger than the p{T)/e{T) ratio accept near the transition temperature, similar as in Ref.[43, 60]. The 
minimum value oi p{T)/e{T) ratio is about .033 in PQMVT model which being slightly larger than the PQM model 
value .026, is less than the lattice result .075 [16, 71]. Similar to the findings of Rcf. [60], interestingly the C^ value is 
found to be less than 0.1 around half the crossover transition temperature in our PQM and PQMVT model results. In 
contrast, using a model of confinement, the results of Ref. [72] find values of about C^ ^ 0.2 around half the transition 
temperature and C^ = 0.15 around the transition temperature. 
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FIG. 10: The variation oi p{T)/e{T) has been shown by thin lines while thick lines show the variation of Cg 



D. Quark number density and Susceptibility 



The first derivative of grand potential with respect to chemical potential gives the quark number density 

_ <90mf 
dfi 



(35) 



The implicit variation of a, $ and $* fields with respect to chemical potential has been accounted for, in the 
evaluation of first derivative of fl{a), Uiog and ft]:^ with respect to chemical potential as given in the appendix A. 
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FIG. 11: Temperature variation of quark number density divided by T^ 



The temperature variation of the quark number density normahzed by T^ in the PQMVT model calculation has been 
shown in Fig. 11 for three quark chemical potentials /i = 280, 294.7 (CEP) and 300 MeV. The dash dotted line shows 
the number density variation for a crossover transition at /x = 280 MeV, here we see a small peak structure. The 
dotted line number density variation shows a sharper rise with a narrow peak at /icEP = 294.7 MeV and it approaches 
the SB value of number density variation (shown by thick dots) for higher temperatures. At fj, = 300 MeV, the Phase 
transition becomes first order, hence the quark number density being a first derivative of the grand potential with 
respect to the chemical potential, shows a jump in the solid line temperature variation. 
The expression of quark number susceptibility is obtained as 



Xq 



d^n 



MF 



dfj,'^ 



(36) 



The second partial derivatives of a, $ and $* fields with respect to chemical potential contribute in the double 
derivatives of ri(cr), Uiog and ftl^ with respect to chemical potential as given in the appendix A. Fig. 12 shows the 
variation of quark number susceptibility normalized by T^ as a function of temperature in the PQMVT model 
calculation for chemical potentials /i = 280, 294.7 (CEP) and 300 MeV. The dash dotted line susceptibility variation 
at fj, = 280 MeV shows a continuous peak structure at the crossover transition temperature. Since at ficEP — 294.7 
MeV, the phase transition turns second order, the dotted line of quark number susceptibility variation shows a very 
large and strongly divergent peak at Tcep. The solid line shows the quark number susceptibility at fi = 300 MeV for 
the first order transition case, we get a discontinuous variation because order parameter registers a jump in the first 
order transition. 



IV. SUMMARY AND DISCUSSION 



We have investigated the temperature variation of the order parameters cr, $, $*, their temperature derivatives 
and various thermodynamic physical observables at non zero and zero quark chemical potentials in the presence 
of renormalized fermionic vacuum term in the effective potential of the PQM model. The results termed as the 
PQMVT model calculations have been compared with the results of PQM model without vacuum term. We have 
used logarithmic Polyakov loop potential for our calculation. 

The chiral crossover transition for the realistic case of explicit chiral symmetry breaking, becomes quite soft and 
smooth at /i = in PQMVT model due to the proper accounting of the fermionic vacuum term contribution in the 
PQM model because the corresponding phase transition at /x = turns second order in the chiral limit of massless 
quarks. The a derivative shows a broad double peak structure at /i = 0. The second peak position at higher transition 
temperature T^ = 186.5 MeV identifies the chiral crossover because the first peak results due to a sharp peak in the 
Polyakov loop temperature variation which signals a rapid confinement-deconfinement crossover transition at /x = 0. 
In a large range of ^, T values (from ^ = and T = 186.5 MeV to /i = 294 MeV and T = 84 MeV), the chiral 
transition remains a crossover and it keeps on becoming sharper with the increase in chemical potential till the point of 
second order transition at hcep is reached in the PQMVT model. Instead of logarithmic form, if we take polynomial 
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form for Polyakov loop potential in our PQMVT model calculation, the temperature derivatives of Polyakov loop field 
<t> and its conjugate $* show distinct non coincident double peak structure in the chemical potential range fi ~ 200 
MeV to ^CEP = 293 MeV and we do not find any double peak structure near /z = in the temperature derivative 
of a field. Hence confinement-deconfinement crossover transition lines for $* and $ fields will be located above the 
chiral crossover phase boundary from /i = 200 to /i = ficEP = 293 MeV. This finding support the quarkyonic phase 
like scenario, having a region of confinement with restored chiral symmetry. Our calculation with logarithmic form 
for the Polyakov loop potential does not support this quarkyonic phase like scenario. Since the chiral transition in the 
massless quark limit is first order at zero chemical potential, the corresponding crossover transition for the realistic 
case has been found to be quite sharp and rapid in the PQM model without any vacuum term. Further the chiral 
transition remains a crossover in quite a small range only from /x = and T^ = 171.5 MeV to /i = 81 McV and 
Tx = 167 MeV in the PQM model calculations. 

The phase diagram together with the location of critical end point (CEP) has been obtained in /i, and T plane 
for Too- = 500 MeV in both the models PQMVT and PQM. The structure of the phase diagram is very sensitive to 
the value of sigma meson mass. For the value ma.= 600 MeV, the transition becomes a crossover in the entire fi 
and T plane for the PQMVT model calculation. We do not have a coincident chiral and confinement-deconfinement 
crossover transitions in the PQMVT model as the chiral crossover transition line lies above the crossover line for the 
confinement-deconfinement transition. Our results of the PQMVT model calculation, are in tune with the standard 
scenario where chiral symmetry restoration occurs at a higher critical temperature than the confinement-deconfinement 
transition temperature. The quarkyonic phase scenario is ruled out in our PQMVT model calculation with constant 
To when we have taken a logarithmic form for the Polyakov loop potential while the calculation with polynomial form 
of Polyakov loop potential supports its occurrence in certain range of /i and T. The critical end point (CEP) gets 
shifted close to the chemical potential axis {^j-cep = 294.7 MeV, Tcep = 84.0 MeV ) in PQMVT model because 
the chiral crossover transition at /i = emerging from a second order phase transition in the chiral limit, becomes 
quite soft and smooth due to the effect of fermionic vacuum contribution in the effective potential and further it 
remains a crossover for large values of the chemical potential. The chiral and confinement-deconfinement crossover 
transition lines are coincident and the CEP of the PQM model gets located near the temperature axis at fJ-cEP = 81 
and Tcep = 167 because the chiral crossover at ^ = 0, having the background of a first order phase transition in the 
chiral limit, is quite rapid and sharp and soon it gets converted to a first order phase transition as we increase the 
chemical potential. 

The temperature variation of thermodynamic observables namely pressure, energy density, entropy density at three 
different chemical potentials (zero, ficEP and fj. > ficEp) has been shown in PQMVT model. Due to the proper 
accounting of appropriately renormalized fermionic vacuum fluctuations, the pressure , entropy density and energy 
density variations at ^ = turn out to be a smoother function of temperature in PQMVT model when it is compared 
with corresponding curves in PQM model calculation. The temperature variations of the interaction measure, speed 
of sound, p{T)/e{T) and specific heat, have been calculated in PQMVT model and QMVT (Quark Meson model 
with vacuum term) model and these results have been compared with the corresponding results in the PQM and 
QM model calculations. Again we find that the presence of fermionic vacuum contribution in effective potential leads 
to the smoother variation of the thermodynamic quantities. Finally we have shown the results of the temperature 
variations of baryon number density and quark number susceptibility at different chemical potentials in PQMVT 
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model calculations. 
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Appendix A: First and second partial derivatives of grand potential 

First partial derivative of logarithmic Polyakov loop potential with respect to chemical potential and temperature 
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Second partial derivative of logarithmic Polyakov loop potential with respect to chemical potential and temperatures 
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Second partial derivative of JlT^ with respect to chemical potential and temperature 



_<^ = -12T 

dy.'^ 



d3p 



(27r) 



s^q^/j 



(A12) 



(A13) 



9T2 



-12 



d^p 



(2^)^ 



2 (B+ + B- 



T D 



q,T 



(A14) 



18 



where 



D„ 



S{B+'+Bf 



9q ' 
9q 



d 



2rT>* 



92$ 



i9x 
92$* 



(-/3Eq 



-;3Eq 



9x2 



9x 
92$ 

9x2 Qy^2 



-/3Eq 



€=-"•* + ^t^=-^°^'1 J^(-^E,*) + ,4J^(-^E,*) 



9x 



9$* 
9x 



-/3Eq 



9x 



(9x 



92x' 



,4|ic-2^E^^^ ^(_^E,-)+ 4-^(-/?E,-) 



(A15) 



C+ = ^e-^^'' + 4$*e-2/'E,+ + 3e-3/3E, 



C- = $*c-'3E." +4$c-2'3E,- +3e-3/^E, 



(A16) 
(A17) 



References 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 
[10] 

[11] 

[12] 
[13] 
[14] 
[15] 
[16] 
[17] 
[18] 

[19] 
[20] 
[21] 
[22] 
[23] 
[24] 
[25] 
[26] 
[27] 
[28] 
[29] 
[30] 
[31] 
[32] 
[33] 



D. H. Rischke, Prog. Part. Nucl. Phys.52, 197 (2004). 
H.Meyer-Ortmanns Rev. Mod. Phys. 68, 473 (1996). 
B.MuUer, Rep. Prog. Phys. 58, 611 (1995). 
L.D.McLerran, B.Svetitsky, Phys. Rev. D24, 450 (1981); 



B.Svetitsky, Phys. Rep. 132, 1 (1986). 



A. M. Polyakov, Phys. Lett. B 72, 477 (1978). 

R. D. Pisarski, Phys. Rev. D 62 111501(R) (2000). 

B. Layek, A. P. Mishra, A. M. Srivastava and V. K. Tiwari, Phys. Rev. D 73 103514 (2006). 
F. Karsch, Lect. Notes Phys. 583, 209 (2002). 

Z. Fodor, S. D. Katz, and K. K. Szabo, Phys. Lett. B 568, 73 (2003). 
A.Ah Khan et al. Phys. Rev. D 64, 074510 (2001). 

C. R. AUton, M. Dormg, S. Ejhi, S. J. Hands, O. Kaczmarek, F. Karsch, E Laermann and K. Redhch, Phys. Rev. D 71, 
054508 (2005). 

Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, Phys. Lett. B 643, 46 (2006). 
F. Karsch, J. Phys. G 31, S633 (2005). 

F. Karsch, e-Print: arXiv:0701.210 [hep-ph]. 

M. Cheng et al, Phys. Rev. D 74, 054507 (2006). 

M. Cheng et al, Phys. Rev. D 77, 014511 (2008). 

S. Digal, E. Laermann and H. Satz, Eur. Phys. J. C 18, 583 (2001). 

J. T. Lenaghan, D. H. Rischke and J. Schaflner-BieUch, Phys. Rev D 62, 085008 (2000). J. T. Lenaghan, D. H. Rischke, 

J. Phys. G 26, 431 (2000). 

S. Chiku and T. Hatsuda, Phys. Rev. D 58, 076001 (1998). 

S. Chiku Prog. Theor. Phys. 104, 1129 (2000). 

T. Herpay, A. Patkos, Zs. Szep and P. Szepfalusy, Phys. Rev. D 71, 125017 (2005). 

T. Herpay and Zs. Szep, Phys. Rev. D 74, 025008 (2006). 

P. Kovacs and Zs. Szep, Phys. Rev. D 75, 025015 (2007). 

G. Fejos, A. Patkos, Phys. Rev. D 82, 045011 (2010). 

O. Scavenius, A. Mocsy, L N. Mishustin, D. H. Rischke, Phys. Rev. C 64, 045202 (2001). 

B. J. Schaefer and M. Wagner, Phys. Rev. D 79, 014018 (2009). 

J. O. Andersen, R. Khan and L. T. KyUingstad, AIP Conf. Proc. 1343 :504-506 (2011) e-Print: arXiv:1102.2779 [hep-ph] 

A. Jakovac, A. Patkos, Z. Szep, and P. Szepfalusy, Phys. Lett. B 582, 179 (2004). 

A. Mocsy, L N. Mishustin, and P. J. Elhs, Phys. Rev. C 70, 015204 (2004). 

B.-J. Schaefer and J. Wambach, Nucl. Phys. A 757, 479 (2005). 

B.-J. Schaefer and J. Wambach, Phys. Rev. D 75, 085015 (2007). 

P. Kovacs and Zs. Szep, Phys. Rev. D 75, 025015 (2007). 

E. S. Bowman and J. L Kapusta, Phys. Rev. C 79, 015202 (2009); J. I. Kapusta, and E. S. Bowman, Nucl. Phys. A 830, 

721C (2009). 



19 



[34] 


A. 


[35] 


L. 


[36] 


P. 


[37] 


J.- 


[38] 


T. 


[39] 


D. 


[40] 


C. 


[41] 


S. 


[42] 


S. 


[43] 


S. 


[44] 


c. 


[45] 


T. 


[46] 


H. 


[47] 


M. 


[48] 


W 


[49] 


K. 


[50] 


K. 


[51] 


K. 


[52] 


H. 


[53] 


P. 


[54] 


K. 


[55] 


A. 


[56] 


K. 


[57] 


B. 


[58] 


B. 


[59] 


B. 


[60] 


B. 


[61] 


H. 


[62] 


U. 


[63] 


G. 


[64] 


V. 


[65] 


T. 


[66] 


R. 


[67] 


A. 


[68] 


M. 


[69] 


J. 


[70] 


L. 


[71] 


S. 


[72] 


B. 



Jakovac and Zs. Szep, Phys. Rev. D 82, 125038, (2010). 
Ferroni, V. Koch, and M. B. Pinto, Phys. Rev. C 82, 055205 (2010). 
Costa, M. C. Ruivo, C. A. de Sousa Phys. Rev. D 77, 096001 (2008). 
L Knuer, M. B. Pinto and R.O. Ramos, Phys. Rev. C 81, 065205 (2010). 

Kahara and K. Tuominen, Phys. Rev. D 78, 034015 (2008); ibid D 80, 114022 (2009). ibid D 82, 114026 (2010). 

Nickel, Phys.Rev. D 80, 074025 (2009). 

Ratti, M. A. Thaler and W. Weise, Phys. Rev. D 73, 014019 (2006). 

Rofiner, C. Ratti, and W. Weise, Phys. Rev. D 75, 034007 (2007). 

Rofiner, T. Hell, C. Ratti, and W. Weise, Nucl. Phys. A 814, 118 (2008). 

K. Ghosh, T. K. Mukherjee, M. G. Mustafa and R. Ray, Phys. Rev. D 73, 114007 (2006). 

Sasaki, B. Friman and K. Redlich, Phys. Rev D 75, 074013 (2007). 

Hell, S. Rofiner, M. Cristoforetti and W. Weise, Phys. Rev D 79, 014022 (2009). 

Abuki, R. Angiani, R. Gatto, G. NarduUi and M. Ruggieri, Phys. Rev D 78, 034034 (2008). 

Ciminale, R. Gatto, N. D. Ippolito, G. NarduUi and M. Ruggieri, Phys. Rev D 77, 054023 (2008). 
.-J. Fu, Z. Zhang and Y.-X. Liu, Phys. Rev D 77, 014006 (2008). 

Fukushima, Phys. Rev D 77, 114028 (2008). 

Fukushima, Phys. Rev D 78, 114019 (2008). 

Fukushima, Phys. Rev D 79, 074015 (2009). 

Hansen, W. M. Alberico, A. Beraudo, A Mohnari, M. Nardi and C. Ratti Phys. Rev. D 75, 065004 (2007). 
Costa, M. C. Ruivo, C. A. de Sousa, H. Hansen and W. M. Alberico Phys. Rev. D 79, 116003 (2009). 

Kashiwa, H. Kouno, and M. Matsuzaki, and M. Yahiro Phys. Lett. B 662, 26 (2008). 

E. Radzhabov, D. Blaschke, M. Buballa, and M. K. Volkov, e-Print: arXiv: 1012.0664 [hep-ph]. 

Fukushima, Phys. Lett. B 591, 277 (2004). 

J. Schaefer, J. M. Pawlowski and J. Wambach, Phys. Rev. D 76, 074023 (2007). 

J. Schaefer and M. Wagner, arXiv:0812.2855 [hep-ph]. 

J. Schaefer, M. Wagner and J. Wambach, e-Print: arXiv:0909.0289 [hep-ph] 

J. Schaefer, M. Wagner and J. Wambach, Phys. Rev. D 81, 074013 (2010). e-Print: arXiv:0910.5628 [hep-ph] 

Mao, J. Jin and M. Huang, J. Phys. G 37, 035001 (2010). 

S. Gupta and V. K.Tiwari, Phys. Rev. D 81, 054019 (2010). 

Marko and Zs. Szep, Phys. Rev.D 82, 065021 (2010). 

Skokov, B. Friman, E. Nakano, K. Redlich, and B.-J. Schaefer, Phys. Rev. D 82, 034029 (2010). 

K. Herbst, J. M. Pawlowski, and B.-J. Schaefer. Phys. Lett. B 696, 58 (2011). 

D. Pisarski and F. Wilczek Phys. Rev. D 29, 338 (1984). 

J. Mizher, M. N. Chernodub and E. S. Fraga, Phys. Rev. D 82, 105016 (2010). 

Quiros, arXiv:hep-ph/9901312. 
Braun, L. M. Haas, F. Marhauser and J. M. Pawlowski, Phys. Rev. Lett. 106, (2011); arXiv:0908.0008 [hep-ph]. 

McLerran and R. D. Pisarski, Nucl. Phys. A796, 83 (2007). 

Ejiri, F. Karsch, E. Laermann, and C.Schmidt, Phys. Rev. D 73, 054506 (2006). 

Mohanty, J.-E. Alam, Phys. Rev. C 68, 064903 (2003). 



